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Abstract 



We analyse the behaviour of the Dirac equation in d = 1 + 1 with Lorentz scalar 
potential. As the system is known to provide a physical realization of supersymmetric 
. quantum mechanics, we take advantage of the factorization method in order to enlarge 

CNj ' the restricted class of solvable problems. To be precise, it suffices to integrate a Ricatti 

equation to construct one-parameter families of solvable potentials. To illustrate the 
procedure in a simple but relevant context, we resort to a model which has proved 
useful in showing the phenomenon of fermion number fractionalization. 



1 Introduction. 



When solving the Dirac equation in d = 1 + 1 with Lorentz scalar potential, the underlying 
supersymmetric structure is crucial. As the system provides a physical realization of 
supersymmetric quantum mechanics (susy qm henceforth), the problem reduces itself to 
a pair of Schrodinger-like hamiltonians related by means of supersymmetry. In doing 
so, both operators share identical energy spectra unless the Dirac equation allows the 
existence of a zero- mode (i.e. a normalizable solution with null energy). Accordingly, we 
would like to show how the methods used in susy qm to obtain analytical solutions of 
the Schrodinger equation can be extended to the Dirac case. As a relevant example of 
the aforementioned models one can consider the problem of the Dirac hamiltonian with a 
scalar potential 4>{x) which characterizes a position dependent band gap. In such a case the 
Dirac equation serves to describe the behaviour of a x-dependent valence and conduction- 
band edge of semiconductors in the Brillouin zone 0. On the other hand, the modified 
Korteweg-de Vries (mKdV) equation exhibits a striking relation to the Dirac equation with 
a time-dependent scalar potential. If (ft(x, t) represents a solution of the mKdV equation, 
the spectrum of the (time-dependent) Dirac operator with the scalar potential <fi(x,t) is 
time independent. In other words, the Dirac operators at different times are unitarily 
equivalent if the potential itself evolves according to the mKdV equation M . 
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Particular interest acquire however the models in which the Lorentz potential can be 
looked upon as the kink-like static solution derived from a classical field theory. When 
going to the quantum level, such solutions give rise in general to new sectors beyond those 
that are seen in the standard perturbative regime. These classical solutions have received 
the names of solitary waves, energy lumps, kinks or solitons. In addition, they appear in 
any number of dimensions: the c/> 4 kink for d = 1 + 1, the Nielsen-Olesen vortex if d = 2 + 1 
or the celebrated 't-Hooft-Polyakov monopole in d = 3+1. (The brave reader can enjoy the 
superb description of the subject contained in the book of Rajaraman ||). Though more 
massive than the elementary excitations, the soliton-like configurations become stable 
since an infinite energy barrier separates them from the ordinary sector. The stability 
is reinforced by the existence of a topological conserved charge which does not arise by 
Noether's theorem from a well-behaved symmetry of the lagrangian, but takes into account 
the large distance behaviour of the field configuration (hence the topological character of 
the solutions themselves). 

On the other hand, dramatic effects appear when a Fermi field is coupled to the 
soliton. Among other things, such models have proved quite useful in the context of the 
fermion number fractionalization which has been observed in certain class of polymers like 
polyacetylene || . Describing the soliton by a non-dynamical c-number field, the behaviour 
of the fermions is governed in this first quantized approach by a Dirac operator which 
includes the background provided by the scalar field. If one considers interactions enjoying 
charge conjugation symmetry, the states of positive and negative energies match between 
them. In addition, the spectrum can exhibit also self-conjugate zero-energy solutions which 
as a last resort assume the responsibility for the fractionalization of the fermion number. 
These results are nicely described in topological terms: no matter what the details of the 
model, the existence of zero-modes is rigorously predicted by general arguments concerning 
differential operators. It may be worth spelling that susy qm represents in its own right 
an alternative tool in the analysis of this connection between physics and topology (see 
further on). 

The organization of the article is as follows. Section 2 takes care of the Dirac equation 
in d = 1 + 1 with static Lorentz scalar potential. There we carefully describe the way in 
which the Ricatti equation enables us to construct one-parameter families of isospectral 
Lorentz scalar potentials for the Dirac case. To illustrate the procedure in a simple but 
relevant context, we resort in Section 3 to the problem consisting in the analysis of the 
Dirac equation on the background provided by the kink-like configuration of the </> 4 model. 
(The first appendix includes a capsule introduction to susy qm while the second one 
outlines the properties of the eigenfunctions associated with the Posch- Teller potential). 

2 Supersymmetric structure of the Dirac equation. 

Although there are several applications of susy qm in the context of the Dirac equation, we 
discuss only the underlying supersymmetric structure of such equation in d = 1 + 1 with 
static Lorentz scalar potential represented by f(x). In doing so, we return to the ideas 
exposed by Cooper et al. some years ago 0. The purpose of this section is to show how 
we can take advantage of susy qm in order to enlarge in a systematic way the restricted 
class of solvable models. To be precise, it is the general solution of a Ricatti equation that 
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provides us with a one-parameter family of potentials which serve to extend the first results 
derived from f(x). To begin with, let us give some words about the notation used in the 
next two sections. Space-time coordinates are represented by x^ (fi = 0, 1; x° = t, x 1 = x). 
As regards the metric tensor we have g^ v = g 1 *", where goo = —gu = 1 and g„v = 
otherwise. Moreover stands for space-time derivatives d/dx^. When considering a 
fermion field ip(x,t) moving in the static Lorentz scalar potential f(x), the behaviour of 
the system is governed by the Dirac equation 

[iVfy - /(z)]#M) = 0. (l) 

Letting 

ip(x,t) = ip(x) exp(— iut) 

the proper Dirac equation reduces to 

j°oj^(x) + i^tp^x) — f(x)ip(x) = 0, 

where the prime denotes as usual derivative with respect to the spatial coordinate. Upon 
choosing now 

we get the pair of coupled equations 

Aip-(x) = ujip + (x), A f tp + (x) = uifj-(x), 

where 

It remains to decouple the equations, i.e. 

A^Aip-(x) = u 2 ip-(x), AA^ip + (x) = u 2 ip+(x) 

which are of the form of susy qm (see appendix A for a brief introduction to the subject), 
namely 

- + [f\x) - f'(x)] rP-(x) = u 2 1>-(x), (2) 

- + [f 2 {x) + /'(*)] = u 2 ^ + (x), (3) 

where the Lorentz scalar potential f(x) is just the superpotential function. To sum up, 
we substitute the Dirac equation on behalf of the two Schrodinger-like equations of a 
genuine susy qm model. In such a case, the components i/>—(x) and ip+(x) represent the 
eigenfunctions of the hamiltonians H- = A^A and H + = AA^ . It is customarily assumed 
that whenever the Schrodinger equations for potentials V± = f 2 (x) ± f'{x) are solvable, 
then there always exists a static Lorentz scalar potential f(x) for which the corresponding 
Dirac equation is also exactly solvable. As expected, the magic of supersymmetry implies 
that both hamiltonians H— and H + become solvable simultaneously since the action of the 
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supercharges establish the right connection between the eigenfunctions ip-(x) and ip+(x). 
The spectrum of energies is of course degenerate except for the hypothetical zero-mode. 

At this point we can recall the way in which the topological behaviour of the superpo- 
tential f(x) dictates the existence (absence) of zero-modes. It is customarily assumed that 
super symmetry, a new kind of symmetry which puts together both bosonic and fermionic 
degrees of freedom, represents one of the most ambitious attempts recently made to achieve 
a trustworthy description of nature. For better or worse, we have no experimental evi- 
dence for such a novel symmetry being exactly realized in physics so that one of the most 
pursued goals in the last two decades has been to describe a suitable mechanism by means 
of which supersymmetry shows up itself spontaneously broken. As susy qm constitutes 
an excellent laboratory to analyse the hypothetical spontaneous breaking phenomenon, 
Witten introduced an order parameter (henceforth the Witten index A) which in general 
serves to shed light on this subject. As a matter of fact the aforementioned index A 
reads [|| 

A(/3) = Tr [exp(-/?ff_) - exp(-/Ji2+)] , 

where (3 stands for a regularization parameter necessary when considering susy qm with 
scattering states. In doing so, A(/3) should represent a measurement of the difference of 
ip-{x) and ip+(x) modes all of them with zero-energy since for positive energies duplica- 
tion occurs between the two parts of the model. In order not to clutter the article, we 
simply point out that when taking into account the subtleties associated with the count- 
ing of states in the continuous part of the spectrum the Witten index A(/3) satisfies the 
equation [|| 

where f± represent the asymptotic values of the superpotential f(x) as x — > ±00. Although 
sometimes cumbersome final indices are obtained if any of the two limits /_ or /+ becomes 
null, it is the case that whenever we observe a change of sign when going from /_ to /+, 
the integration of (||) gives 



A(/3) 



being $ the Fresnel probability function. Now it suffices the limit j3 — > 00 to reveal the 
topological roots of the model: the Witten index amounts to ±1 so that the asymptotic 
behaviour of f{x) dictates the existence (absence) of the non-degenerate zero-mode (see 
Appendix A). From a qualitative point of view, this result can also be seen as follows: a 
normalizable zero-energy eigenfunction exists for the superpotentials just exhibiting the 
right behaviour at infinity. In summary, the odd-like superpotentials yields a zero- mode 
while the even-like ones are not capable of producing such a state. 

However, the above discussion is only part of the whole story. Once we assume 
the existence of a solvable quantum mechanical problem for the pair of one-dimensional 
Schrodinger-like equations written in @) and (§), different procedures can be profitably 
used for generating one-parameter families of models which share exactly the same set 
of eigenvalues. The reasoning becomes much more accesible when resorting to susy qm 
in the spirit of the factorization schemes exposed by Mielnik @. To fix the ideas, let us 
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take /+ > and /_ < so that a well-behaved zero-mode exists associated with H-. In 
doing so, the energy spectrum of H + is positive whereas starts from E = 0. Being 
tjj-o(x) precisely the zero-energy eigenfunction of H_, the static Lorentz scalar potential 
f(x) appears as the logarithmic derivative of such a state, i.e. f(x) = —iP'-o(x)/i/j-q(x). 
At this stage, we ask ourselves whether the factorization H+ = AA^ is unique or can be 
more general than usually realized. Writing now a new factorization like 



for 



B 



>l +F(x), = —^ + F(x) 



dx dx 
the function F{x) must satisfy the Ricatti-like equation 

F'{x) +F 2 (x) = V+(x). 



(5) 



Whenever we know a particular solution of (|5|) (the potential f(x) itself!), the general 
solution can be obtained by putting || 



F(x) = f(x) + <f>(x) 
so that 

<t>'{x) +<p 2 (x) + 2f(x)<f>(x) 



0. 



Now it proves convenient the introduction of p(x) 
first-order differential equation like 

-p'{x) + 2f{x)p{x) + 1 = 

whose general solution corresponds to 



(6) 

l/4>(x) to transform @ into a 



p(x) 



exp 



2f(t) dt 



A + 



exp 



2f(z) dz dt 



where in principle A is a real number. Bearing in mind the relation between f(x) and the 
zero-mode if)'-o(x), we get for the function F{x) the useful expression 



F(x) = f{x) + 



^o(^) 



whenever the choice for A avoids any singularity along the real axis. To sum up, we go 
from the small potential f(x) to the large potential represented by F{x) just by solving 
in the standard way a Ricatti differential equation. In doing so, we can take advantage 
again of the commutation formula of susy qm to get a family of hamiltonians H , i.e. 

H = B^B 

whose close relation with if_ can be explained as follows. Except for the zero-energy 
eigenstate, the potential V(x) included in H, namely 



V(x) = V+(x) 



I d 

dx 
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becomes isospectral with respect to the V + {x) itself. As anticipated in the second section, 
the zero-modes are found from first-order differential operators so that in our case it 
suffices to consider 



whenever A < — 1 and A > in order to guarantee the right normalization of the state. 
In this regards, the procedure may be understood as a renormalization of the zero-mode, 
leading to a one-parameter of hamiltonians H isospectral with the original i2_ . Returning 
to the Dirac equation itself, the situation is both surprising and unexpectedly subtle. 
First of all, we start from a static Lorentz scalar potential f(x) which yields solvable 
Schrodinger-like potentials V±(x) = f 2 (x) ± f'(x). Afterwards, we enlarge in a systematic 
way the restricted class of such solvable models just by going from the small potential f(x) 
to the large potential F(x) through the right combination of the commutation formula of 
susy qm and the general theory of Ricatti equations. In addition, the down components of 
the spinor remain unchanged whereas the transformation ip-(x) — > ip{x) must be carried 
out as far as the upper components are concerned. In any case, we avoid the difficulties 
associated with the handling of the potential V(x) since the ip{x) themselves can be 
obtained from ip+(x) by the action of the operator . 

3 Illustration: kinks in the presence of Fermi fields. 

To illustrate the above formalism in a simple but physically relevant scenario, we resort 
now to the Dirac equation in the background provided by the kink-like configuration of 
the <j) A model. Such system has proved quite useful in the context of the so-called fraction- 
alization of the fermion number which can be observed in polymers like the polyacetylene. 
By any standard, the search for classical solutions represents nowadays one of the most 
prolific areas in quantum field theory. The interest is mainly motivated by the belief 
that such classical configurations, which of course make stationary the action, may shed 
light on the properties of the system at issue. In this section we restrict ourselves to 
non-dissipative configurations with finite energy because they constitute the natural can- 
didates to describe new sectors of the model beyond the perturbative regime. On the 
other hand, dramatic effects like the fractionalization of the fermion number appear when 
a Dirac field is introduced to the kink itself. At this stage, it should be interesting to re- 
view the pioneering model of Jackiw and Rebbi || in the light of the formalism exposed in 
the previous section. For reference, we recall that the behaviour of the system is governed 
by the Lagrangian density 



where as usual 4>(x) and ip(x) stand for the boson and the fermion fields respectively. 
Notice that for the sake of simplicity the coupling constants have been absorbed in the 




to obtain 




C = ~ d^d^ - l - (</> 2 - l) 2 + - 



(7) 
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fields themselves. Unless otherwise noted, it suffices here to discuss a first quantized 
approach where the fermion moves in the non-dynamical c-number field provided by the 
static kink 4>k{x). In doing so, the Dirac equation derived from (|7|) reads 

i^d^(x, t) - 2(j) k (x)^{x,t) = 

which is of the form written in (|l]). In other words, the kink-like profile represents in 
this case the static Lorentz scalar potential. At this stage, it is customarily assumed the 
existence of the well-known kink of the 4> 4 model, namely || <f)k{x) = tanhx a topological 
configuration which smoothly interpolates between the constant configurations 0_ = — 1 
and <p + = 1. On comparing now with the general formalism of the second section, we 
face a standard Dirac problem with f(x) = 2 tanh x whose properties (spectrum and 
eigenfunctions) are well studied (see Appendix B). For reference, we simply point out that 
in this case we have 

F_(x) = 4 V+(x) 



cosh 2 x ' cosh 2 x 

while the zero- mode if)-o(x) reads 

y/3 1 



4>-0{x) ~ „ , 2 ■ 

l cosh x 
Now we can write that 

F(x) = 2 tanh x + 



3 cosh 4 . 



4A + 3 tanh x - tanh 3 x + 2 ' 

As regards the one-parameter family of potentials V{x) we have 

~ . . 6 d 3 cosh -4 x 

V{x) =4 ^ 2— 

cosh x dx 



4A + 3 tanh x - tanh 3 x + 2 



while the renormalization of the zero-mode translates into 
2^/3 a/A(A + 1) cosh -2 x 



ip (x) 



4A + 3 tanh x - tanh 3 x + 2 



The procedure closes itself when writing the spinor solutions of the large potential F(x) 
in terms of the components if)+(x) and ip{x) according to the general method. 



4 Conclusions. 

In this article we have considered in detail the procedure for obtaining exact energy eigen- 
values and eigenfunctions for the Dirac equation involving a Lorentz scalar potential in 
d = 1 + 1. The generic model is useful in different branches of physics, ranging from 
condensed matter to nonlinear equations. For instance, it serves to study the fermion 
number fractionalization observed in certain class of polymers like polyacetylene, the po- 
sition dependent band gap in semiconductors or the theory of solitons of the modified 
Korteweg-de Vries equation. Taking into account the way in which the system itself pro- 
vides us with a realization of super symmetric quantum mechanics, we face ultimately a 
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pair of Schrodinger-like hamiltonians. As expected, the two aforementioned operators 
share the energy spectra except for the hypothetical existence of a zero-mode, i.e. a 
normalizable solution with E = 0. 

In addition, the so-called topological techniques have proved quite useful in modern 
physics. When we need to analyze the spectrum of a well-behaved linear differential 
operator, the zero-modes are the most crucial ones for applications. To determine the 
properties of these modes we can resort to the results contained in the index theorems: no 
matter what the local details of the model at issue, the existence of zero-modes depends 
on the global behaviour of the system (hence the topological character of the procedure 
itself). Taking advantage of the structure associated with susy qm, it suffices to consider 
a Ricatti equation in order to enlarge the restricted class of problems whose eigenvalues 
and eigenfunctions can be presented in closed form. To sum up, whenever we have a 
pair of solvable Schrodinger-like equations then there also exists a one-parameter family 
of Lorentz scalar potentials for which the Dirac equation is exactly solvable. 



Appendix A 

In this appendix we present a brief introduction of one-dimensional susy qm, as tailored 
to our needs. In its most simple formulation the model deals with a time-independent 
hamiltonian matrix H s of the form 



H- 
H+ 



which contains both one-particle hamiltonians iT_ and H + . This operator H s is part of a 
graded algebra (superalgebra in the jargon of physics) whose closure is achieved by means 
of the right mixing of commutation and anti-commutation relations. At this point, it is 
customarily assumed the existence of the so-called superpotential function W(x) which 
gives rise to the couple of first-order differential operators A, A^ given by 

A = ^L + W(x), = — ^ + W(x). 

ax ax 

Now we can write the supercharges Q, Q\ i.e. 

which together with the hamiltonian H s constitute the superalgebra of susy qm, namely 

H S = {Q,Q^}, [H s ,Q]=[H s ,Q^=0, {Q,Q} = {Q\Q^} = 0. 

In more physical terms susy qm can be understood as a matrix hamiltonian H s acting 
on two-component wave- functions ip(x) like 

'ip-(xY 



ip(x) 



ip+(x) 



where it should be emphasized that H + is obtained from H- just by reversing the order 
of the operators A and A^ . In other words 

H- = A^A, H+ = AA^ 
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so that the two one-particle hamiltonians H- and H + read 
for 

V±(x) = W 2 (x) ± W'(x), 

where the prime denotes derivative with respect to the spatial coordinate. It is often 
the case to refer to the pair of potentials V-(x) and V+(x) as supersymmetric partners 
in accordance with the standard susy qm language. To finish this brief exposition, let 
us outline the main properties of our model as far as the standard realization of susy 
qm is concerned. First of all, the energy eigenvalues of both H_ and H + are positive 
semi-definite since the eigenkets ip-(x) and ip+{x) verify in fact the relations 

<V-|#-IV+) = <V>+|tf+IV>-} = \\A^-)f, 

where |||'i/'±)|| represents the norm of the eigenstate \ip±}- As regards the hypothetical 
zero-energy solutions ^_o(x) or iP + q(x) we have that 

^_o(x) = 0, Aty+ofc) = 0. 

In doing so, they can be obtained by solving first-order differential equations so that we 
write the formal solutions 

i/> M (x) ~ exp (± f W{y) dy^j . (8) 

As any well respected normalizable zero-energy eigenfunction requires boundary condi- 
tions of the type ip±o — > as x — ► dboo, the meaning of (§) is clear: the zero-mode, if any, 
must be non-degenerate. The above argument illuminates the framework in which super- 
symmetry itself appears either unbroken or spontaneously broken. To be more specific, if 
susy remains exact the energy of the non-degenerate ground-state is zero while the spon- 
taneous breaking phenomenon translates into twofold ground-state. On the other hand, 
all other energy eigenvalues appear duplicated. When considering the eigenket ip- (x) with 
energy i.e. 

the action of the operator A on the t/j- (x) itself provides an eigenfunction of the partner 
hamiltonian H + . In other words 

H+ [Aip_(x)] = AA^Aip-tx) = E_ [Aiff-(x)] . 

Just by reversing the order of the operators, we start now from the eigenkets ip+(x), 
namely 

thus obtaining as expected eigenfunctions of i?_ through the action of A< , i.e. 



H- 



^^+(2;) = A^AA^ + {x) = E + Aty+(x) 
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To sum up, except for the hypothetical zero-modes all other energy eigenvalues appear 
duplicated when the whole model is studied. Further on, are the first-order differential 
operators A, A' that establish the connection between the two sectors of the system. 
However, the above structure exhibits important advantages if we study the problem 
from a practical point of view. If we resort for instance to a computational technique to 
estimate the energy spectrum, the original hamiltonian can be substituted on behalf of 
the isospectral one to improve the properties of convergence. The same argument can be 
applied of course when considering variational methods. 



Appendix B 

For reference, we outline in this appendix the mean features concerning the Posch- Teller 
potentials which arise in the analysis of the Dirac equation on non-trivial backgrounds 
like the one associated with the kink of the (j) 4 model. The problem itself can be elegantly 
solved in a Lie-algebra framework according to the technique described by Frank and 



Wolf [10], whose notation will be used in the following. To start from scratch, let us 
consider the Posch- Teller hamiltonian 



+ 



2 dx 2 cosh 2 x sinh 2 x ' 

where we choose a system of units in which h = m 
introduction of the parameters ki (i = 1, 2) given by 



(9) 

1. Now it proves convenient the 



(2h 



1) 



1 

s = - 
2 



(2k 2 - If 



which in turn enable us to write the energies E^ of the bound states associated with the 
Posch- Teller potential of (^), namely 



where k - 



(2k - 1) 



&min, fcmin-1, • • • > 1/2 with k 



(10) 



k± — k 2 . Letting k to be of the form 



I + IK 
k = k > 



the expression of in (|T(]) also takes care of the continuous part of the spectrum. It 
remains to identify the behaviour of the wave-functions themselves which, as expected, 
depends on the values of the singularity parameter s in @. At first glance the situation 
is as follows: 

i) For the bound states we have 



(ki,k 2 ) 



(coshx) 



-2fci+3/2 



(sinhx) 2fc2 ~ 1/2 F[a,/? : 



r,z\ 



where F [a, /3, 7; z\ stand for the 2-F1 Gauss hypergeometric function with the identifica- 
tions 



a = —k\ + k 2 + k, (3 = —k\ + k 2 — k + 1, 7 = 2k 2 , z = — sinh 2 x. 
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ii) When considering the scattering states we find 
now with 

a = k\ + k 2 — k, (3 = k± + k 2 + k — 1, 7 = 1k 2 , z = — sinh 2 x. 

With regards to singular potentials like the ones written in (|9|), the careful analysis of 
the problem can be carried out in accordance with the subtle mathematical tools exposed 
by Reed and Simon [11]. As far as our model concerns, it might seem plausible to infer 
the existence of three well different regions depending on the values of the parameter s, 
namely: 

i) If s < —1/8, the model itself becomes physically unrealistic since the spectrum is 
not bounded from below. 

ii) Whenever —1/8 < s < 3/8, the singularity is not strong enough to make the wave- 
functions vanish at the origin. However, it proves necessary to define one-parameter 
family of self-adjoint extensions so that a bounded from below spectrum appears 
once a proper domain is specified. In more physical terms, the wave-functions pass 
across the singularity point and the model extends itself along the real axis. 

hi) Finally, if s > 3/8 the singularity acts as an impenetrable barrier thus dividing the 
space into two independent regions, i.e. x > and x < 0, forcing the wave-functions 
to vanish at the origin. 

The above results shed light on the family of susy qm models governed by superpoten- 
tials like 

W(x) = £tanhx, £ = 1,2,... (11) 
which give rise to 

h. = -— + e - H+ = -—+e- 

dx 2 cosh 2 x ' dx 2 cosh 2 x 

In accordance with the previous results we have for H_ 

I 3 , £ 1 
fcl = 2 + 4' ^ 1 = "2 + 4 
together with 

k 2 = ~, k 2 = ~ 

so that the discrete energy spectrum reads 

E-j =£ 2 -(£-j) 2 , j = 0,l,...,£-l 
while the scattering states start at E = £ 2 . With regards to H + we find 
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and again 





with energies 



E. 




>2 




j = !,...,£-!. 



'+. 



The difference between both discrete spectra reduces itself to the zero-mode of iJ_ 
which adopts the form 

ip-o(x) ~ P/(tanhx), 

where the P/(tanhx) represent associated Legendre polynomials. Bearing in mind that 



the logarithmic derivative of the zero-mode itself provides us with the superpotential 
functions written in (|TT|), i.e. 



In this case, being s = 0, we get rid of the difficulties associated with the singular 
behaviour of the potential. Among other things, the system extends along the real line 
and the expected pattern of even and odd wave-functions (notice that we are dealing with 
even potential) appears in the end. 
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Pg(t&nhx) ~ cosh * x 



W(x) 



P'Ktanhx) 
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